Shapes of Semiflexible Polymers in Confined Spaces 
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We investigate the conformations of a semiflexible polymer confined to a square box. Results of 
Monte Carlo simulations show the existence of a shape transition when the persistence length of 
the polymer becomes comparable to the dimensions of box. An order parameter is introduced to 
quantify this behavior. A simple mean-field model is constructed to study the effect of the shape 
transition on the effective persistence length of the polymer. 



INTRODUCTION 



Biological macromolecules live in the crowded, con- 
fined environment of a cell. They are often packed into 
spaces that are much smaller than their natural size and 
adopt conformations which are unlikely to occur in free 
space. For example, viral DNA is packaged in a cap- 
sid whose dimensions are comparable to the persistence 
length of DNA [H, 0, i, i, d, which obliges the DNA to 
adopt a tightly bent shape. Other examples of strong 
confinement include actin filaments in eukaryotic cell Q , 
and protein encapsulated in E.Coli f?\. In such environ- 
ments, the macromolecule is forced to optimize the com- 
peting demands of configurational entropy, excluded vol- 
ume and bending energy. The packaging of DNA in viral 
capsids has been investigated extensively and illustrates 
the effects of these competing interactions and the addi- 
tional electrostatic forces arising from its charged nature 
[li, ^2,] . Experimental and theoretical studies 3., ^, ^] show 
that the static, and dynamic properties such as thresh- 
old force for packaging and packaging time, change dras- 
tically with changes in confining geometry. Another ex- 
ample is provided by driving DNA through a nanochan- 
nel by electrical force [3, [lH , where the dimensions of 
the channel determine the translocation time and thresh- 
old voltage. Understanding the statistics of polymer in 
confinement is therefore essential for exploring both the 



static and dynamic properties of DNA [12 1 



In this paper, we investigate the shape and rigidity of 
a single, semiflexible polymer under confinement. Using 
the Bond-Fluctuation-Model (BFM) [Tsl, we simulate a 
single, semiflexible polymer chain confined to a box in 
two dimensions. One advantage of the two-dimensional 
(2D) simulations is that they can be compared directly to 
experimental measurements, since conformations of DNA 
on a substrate can be visualized through atomic force mi- 
croscopy (AFM) which characterize the statistics of fluc- 
tuations of DNA[i^. Our simulations show that there 
is a shape transition as a function of increasing rigid- 
ity. We analyze the nature of this transition by defining 
an order parameter associated with the transition. As 
expected, the flexibility of the chain, measured by the 
tangent-tangent correlation function, evolves with the 



shape change. We construct a theory for this correlation 
function by considering Gaussian fluctuations around the 
"average shape" . 



SIMULATIONS 

The simulations are performed using the bond fluctu- 
ation algorithm, introduced by Kremer and Carmersin 
[l^ . The bond fluctuation model (BFM) is a coarse- 
grained model of a polymer in which the chain lives on 
a hypercubic lattice and fluctuations on scales smaller 
than the lattice constant are suppressed. The polymer is 
represented by a chain of effective monomers connected 
by effective bonds. The effective bonds and effective 
monomers are constructed to account for excluded vol- 
ume effects. In 2D, each monomer occupies an elemen- 
tary square and forbids other monomers to occupy its 
nearest and next-nearest neighbors . In order to avoid 
crossing of bonds, the effective bonds are restricted to the 
set: 



B = P(2, 0) U P(2, 1) U P(2, 2) U P(3, 0) 
UP(3,1)UP(3,2) 



(1) 



in units of lattice spacing, and P{m, n) stands for the sign 
combinations ±m, ±n and all the permutations. With 
this restriction, the average bond length (Kuhn length 
b) is 2.8. Fig. [T] illustrates the BFM representation of 
a polymer and shows one allowed move that obeys the 
self-avoidance condition. 

The only non-bonded interaction included in this rep- 
resentation is excluded volume. To represent a semiflexi- 
ble polymer, such as DNA in this model, the energy cost 
of bending has to be incorporated. The mechanical prop- 
erties of DNA are well described by the worm-like chain 
model (WLC) , in which DNA is characterized by its 
total contour length L and persistence length Ip. In the 
continuum limit, the bending energy Hi, in D-dimcnsional 
space can be expressed as: 



ds 



(2) 
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FIG. 1: One configuration of polymer in BFM. The bond 
vectors are restricted to the set B (Eq. [l| and one possible 
move is denoted by the dashed line. 



where the stiffness, k, is related to Ip by = ^^^j^lp, 
and Ip characterizes the decay length of the correlation 



between the tangent vectors [l^ 
(u(s) • u(0)) cx e" 



(3) 



Here u(s) = 9R(s)/9s is the tangent vector at arclength 
s, and R(s) is the position vector. The effect of confine- 
ment on the form of this correlation function is one of 
the questions we explore in this paper. In the following, 
all energies will be measured in units of ksT. 

In a lattice representation such as the BFM, the bend- 
ing energy Hi, associated with a configuration is ex- 
pressed as [it} : 



N-l 

26 ^ 

i=l 



where N is the total number of monomers in the chain. 
The modified BFM model now includes excluded volume 
interactions and bending rigidity. Monte Carlo calcu- 
lations are carried out using the Metropolis Algorithm 
with Hh as the energy function since excluded volume 
effects are incorporated in the allowed configurations of 
the BFM. 

Fig. [2]illustrates the simulation setup for a semiflexiblc 
chain in a square box with linear dimension W, in which 
all lengths are measured in units of the Kuhn length. 



u(s) 




W=15 
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FIG. 2; (a) One initial configuration of our simulation. The 
contour length of the chain is 60, box size is 15. u(s) is 
illustrated by arrows. (6) Mapping part of the chain circled 
in (a) onto u(s) space. 



15 and L — 40, = 10 to investigate the effects of L 
and W. In addition, we compare the results of a circular 
box to a square box. Each Monte Carlo trajectory is 
allowed to equilibrate for 10^ Monte Carlo steps which 
correspond to many (~ 100) Rouse relaxation times [l^. 
The range of the persistence length Ip is varied from 0, 
corresponding to a flexible chain, up to W, which is the 
rod- like limit. For each more than 20 independent 



The radius of gyration Rg of the unconfined semiflexible runs were performed to get adequate statistics, 
chain is proportional to L^lp~'^ [l8|, where v is the Flory 
exponent (z^ = | for 2D, v w 0.588 for 3D 19]) . We 

set up conditions such that Rg is larger than the box SHAPES IN CONFINEMENT 
size. For a chain of contour length L — 60, W has to be 

less than 21.5, for this condition to be met. Most of our For a semiflexible chain, the mechanical properties are 

results were obtained for L = 60 and W = 15, however, well described by a self-avoiding walk [23|, and the only 

limited sets of data were also collected for L = 80, W = effect of the semiflexibility is to change the Kuhn length 



3 



|18l |. In confinement, especially if the confining dimen- 
sions are smaller than or comparable to the persistence 
length, semiflexibility plays a much more important role 
since the bending energy is in intense competition with 
the entropy. As a result, the conformations adopt an av- 
er a_ge_shape which is influenced by the confining geometry 
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We visualize and classify a chain's configurations using 
the average tangent vector, (u(s)), along the chain. The 
chain configuration is mapped onto the tangent vector 
space by calculating the tangent vector at each monomer, 
and then connecting the end points of these vectors ac- 
cording to their chain sequence, as shown in Fig. 2(b) 



The average (u(s)) can then be calculated by averaging 
over Monte Carlo trajectories. We choose (u(s)) instead 
of position vector to characterize the shape because the 
former is independent of the starting point of the chain. 
We average over 10^ statistically independent configu- 
rations generated by a Monte Carlo trajectory. Fig. [3] 
shows (u(s)) for = 1 and Ip = 5. These two shapes 
look very similar: tangent vectors are arranged randomly 
so that the connecting bonds cross each other frequently. 
This is not surprising since for Ip << W,L, the confor- 
mations are close to that of a self-avoiding walk and there 
is no "order" in the tangent vector space. 

The situation is dramatically different for /p > as 
illustrated in Fig. [H For 5 < < 15, there are vestiges 
of order in u(s) space but for Ip > 15, the configurations 
demonstrate clear signatures of ordering as measured by 
(u(s)) (Fig. |4]). At the onset of ordering around Ip ~ W, 
there are extended correlations between the tangent vec- 
tors but the shape is more like an ellipse than a circle. We 
believe that the reason behind this is the square shape of 
the box. In the box, the diagonal direction offers maximal 
space, therefore when the bending energy first starts se- 
riously competing with the chain entropy, an anisotropic 
shape can better optimize the entropy than an isotropic 
one. To check the validity of this idea, we analyzed the 
shape of a chain confined in a circle with similar param- 
eters: L = 60, Ip = 15 and a diameter of 15. In this case, 
(u(s)) was found to be isotropic. For Ip much larger than 
the box size there is a clear spiral shape to the chain (Fig. 

EM- 

In order to better quantify the order in the shapes 
and to explore the nature of the transition from a state 
with no visible ordering to the clear spiral shapes for 
Ip > 15, we define an order parameter, 5" associated with 
the tangent vector: 



OS 



(4) 



Since the chains live in 2D, '5 can be considered as a 
scalar order parameter. From geometry [23:|, the direc- 
tion of u(s) x^is along the binormal vector and its 
magnitude is the curvature. Therefore ^' is zero for a 




(a) 





(b) 



FIG. 3: (a) Plot of (u(s)) using the mapping illustrated in 
Fig. |2(b)| for a polymer of Ip — 1. (6) Ip = 5. 



shape with no preferred bending direction such as an un- 
dulating line but is nonzero for a shape such as a spiral. 
Measuring j^*! as a function of Ip for a given L and W, 
we find a sharp change at /p ~ 15 which is close to the 
box size W, as shown in Fig. [5] 

To analyze the nature of this transition, we measured 
the distribution of the order parameter at different points 
along the transition. 

At ?p = 5, the distribution is close to a Gaussian, as 
seen in Fig. [S] With increasing values of Ip the distri- 
bution becomes significantly non-Gaussian (Fig. [7]) and 
ultimately evolves to a bimodal distribution resembling 
two Gaussians(Fig. [8|). The distributions are reminis- 
cent of a first order transition, except for the fact that 
for Ip ~ W, the distributions clearly show four peaks as 
evident in Fig. [T] 

Typical configurations corresponding to the outer 
peaks are similar to the plot in Fig. [9l spiral in a box. 
A typical configuration for an inner peak is illustrated in 



Fig. 10(a) and shows an elliptical character. In a circular 
box, the inner peaks are still evident but they do not have 
the elliptical character (Fig. |10(b)[ ). These peaks are. 



(a) 



FIG. 5: Absolute value of the order parameter, l^*], as a func- 
tion of l„. 
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FIG. 4: (a) Same as Fig. |3]but for Ip = 15. (6) Ip = 25. 



therefore, more robust than the eUiptical shape. Since 
the only other length scale in the problem, besides, Ip and 
W is the contour length or, equivalently, Rg, we have fu- 
ture plans of exploring the distributions as a function of 
Rg in order to better understand this intermediate state. 

Interactions of a semiflexible polymer with hard (or 
semisoft) boundaries are known to introduce orienta- 
tional correlations [U [23, HJ]. One example is the 
transition to liquid crystalline phases of polymers con- 
fined between two hard walls [221 when the persistence 
length is comparable to the dimensions of the confine- 
ment. The shape changes that we are observing in our 
study have a similar physical origin. 



TANGENT- TANGENT CORRELATIONS AND 
EFFECTIVE PERSISTENT LENGTH IN 
CONFINEMENT 

^From the discussion above, it is clear that the shape, 
characterized by the average tangent vector, depends on 



FIG. 6: Distribution function P(*) for Ip = 5. The red line 
denotes a Gaussian fit. The mean is at = and the config- 
urations are disordered (Fig. |3}. 
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FIG. 7: Distribution function P(vE') for Ip = 15. There are 4 
peaks and a valley at '3/ = 0, and this is no longer described 
by a Gaussian. 
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the relation between Ip and W. The tangent-tangent 
correlation function, C(s, s'), defined as 



C{s, ^ ((u(.) - (u(,s))) . (u(.') - (u(.')))) (5) 

is expected to depend on the shape of the polymer. Such 
changes have been observed in experiments. For ex- 
ample, in actin filaments trapped in narrow channels, 
C(s, s') exhibits a oscillatory behavior, and an effective 
persistence length, deduced from the correlation func- 
tion, shows changes from the bare persistence length 
[^, [2^ . Because of translational invariance and the indis- 
tinguishability of two ends of a polymer chain, C(s, s') 
is only a function of |s — s'|, and in the simulations, we 
measure C(s) = C(s,0). Figs, [m and [H] show C(s) 
measured at Zp < 5 and 8 < Ip < 18, respectively. 

For Ip < 5, C(s) exhibits exponential decay with a 
length scale which is indistinguishable from the bare per- 
sistence length, Ip. In this regime, where Ip is much 
smaller than W, the free energy of the polymer is dom- 
inated by the entropy , and the fluctuations resemble 
those of an unconfined semifLexible polymer. 



C(s) 



1:1-5 

p 



"'■•*-•"— I =5 



10 20 30 40 50 



FIG. 11: C{s) (EqEl) for Ip < 5. When Ip < 5, C(s) decays 
exponentially and then fluctuates around 0; a behavior that 
is similar to an unconflned chain. When Ip = 5, C(s) begins 
to show the oscillatory structure. The meansured persistence 
length le ^ Ip. 
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FIG. 12; C(s) (EqlSll for 8 < Ip < 18. Graphs are labelled by 
the value oflp. In contrast to the behavior observed at small 
values of Ip (fTTj). there is a distinct periodic behavior. 



The situation changes for Ip larger than or comparable 
to W. As seen from Fig. [T^ C{s) exhibits an oscillatory 
character. A negative value of C(s) indicates a reflection 
in orientation of the tangent vector. The oscillations, 
therefore, reflect the constraints imposed on C(s) by the 
geometry of the confinement and we expect the periodic- 
ity to be related to the size of the box, W. We therefore 
represent C{s) as Ae~T^ cos -g, with A,B and le as fitting 
parameters. From Fig. [T3land[T4l we see that this form 
captures the properties of C(s). In the regime Ip > 8, 
the parameter B is not sensitive to Ip and approaches 
the value 7.5, which is equal to W/2. We can extract the 
effective persistence length, le, from the fits. The vari- 
ation of this length, which characterizes the rigidity of 
the polymer under confinement is shown in Fig. 1151 It is 
clear that 1^ begins to deviate from Ip for Ip > 8. This is 
the regime in which the distribution of the order param- 
eter starts exhibiting multiple peaks but the average 
order parameter is still close to zero. 



GAUSSIAN THEORY OF C(s) 

To understand the origin of the changes in C(s) with 
the evolving shape of the polymer, we construct a Gaus- 
sian model of the fluctuations of the tangent vector. In 
this model (Fig. [T6|) , the polymer fluctuates around a cir- 
cle with radius d, and the magnitude of the fluctuation 
(j){s) is assumed to be small. To describe the energet- 
ics of the fluctuations 0(s), the bending energy term of 
the Hamiltonian for a worm-like chain is augmented by a 
harmonic conflnement potential Hj, that penalizes devia- 
tions from the circle: Hj = ^ 4''^{s) ds. The harmonic 
potential mimics the confinement effects and the param- 
eter, A, is an effective coupling constant. This type of 
modeling has been successfully applied to analyze a poly- 
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FIG. 13: Fit of C(s) by the function Ae cos -| for Ip = 8. 
Box size ly = 15, and .B ~ -Si^. 
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FIG. 14: Fit of C(s) by the function Ae cos -| for Ip = 25. Ze 
changes and B ~ 




FIG. 15: Plot of the effective persistence length U versus the 
bare persistence length Ip for W = 15, L = 60. In the inset, 
the black dots represent le vs. Ip for Ip > 13 and the red 
triangles represent the plot of the function (I20|l . 
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FIG. 16: Mean field model for calculating tangent-tangent 
correlation function. 



mer confined in a tube 



H = 



fit 



2a |. The total Hamiltonian is: 

L 



du(s)-|2 , A 



4f{s)ds 







Introducing two unit vectors, eg{s) and ep(s) 

ee(s) = {-sin- cos-) 
a a 

ep(s) = {cos-, sin^) 
da 

the tangent vector can be expressed as 

d4){s) 



u{s) = {l + ^c^{s))eg{s) + 



-ep(s) 



(6) 

(7) 

(8) 

(9) 
(10) 



ds -^"^ 

Then the Hamiltonian H takes the following form 
kL f^,, K Ad^ , ,,, , „ ,dib(s).^ 

where "0(5) — the first term represents the bending 
energy of a circle and the integral represents the energy 
from the fluctuations. Here we make use of /J" 'ip{s) ds = 
since the length of the polymer is fixed. In addition, 
we have retained only the leading order gradients of ijj 
in the harmonic potential part in order to be consistent 
with the terms in the WLC model. The Green's function 
is defined as G{s) = (V'(s) • ^-(0)). 

In the limit of L — > cxd , we can decompose the fluctu- 
ations into Fourier modes, 



'0(9) dq 



(11) 



and obtain an expression for H and G{s, s') in the Fourier 
space, respectively. 



G{q)^ 



27r 



(13) 



The Green's in real space is then given by G{s) oc e 
where the length If. is: 



6k 



Ad2 



(14) 



The correlation function G{s) is related to G(s) by 

^2 



\s\ ,2 d 



C(s) = cos^(l-d^-^)G(s) 
d os^ 



-'2d sin 



\s\ dG{s) 



d ds 

Using the form of G{s), we obtain: 



C{s) 



. \s\ 



d' \s\, 



(15) 



(16) 



In the limit of small f- 



the leading contribution to 

C(s) is e~'^cos|, a form that is in remarkably good 
agreement with the simulation results shown in Figs I13[ 
[T4l for which the condition d << le is satisfied. The 
radius d of the circle can be obtained by minimizing 
the mean-field free energy [i^], F, corresponding to the 
Hamiltonian H: 



(17) 



Integrating out the Gaussian functional, leads to: 

Ad2 



F ■ 



kL 
2dP 



L 
2^ 



Hinq' + ^+--^)dq (18) 



Setting |j = 0, yields 



12^2 



Ad4 



(19) 




The effective coupling constant A was introduced to 
mimic the effects of the confining box. Making use of 
the observation that d = in the simulations, we can 
determine A in terms of W and k. Using this result, we 
subsequently obtain If. in terms of Ip and W: 



V3 
2 



M/2 



611 



2V3L 



ill 



(20) 



With increasing values Ip, le is predicted to increase lin- 
early with Ip for Ip ^ Comparing to the results of 

simulations fFigfTS]). it its seen that the predictions of 
the Gaussian theory are in semi quantitative agreement 
with the simulations for Ip > 13. This linear dependence 
is different from the stiff polymer confined in the tube, 

for which le(xll [i,|2|. 
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SUMMARY 

We have investigated the conformational and elastic 
properties of a single semiflexible polymer confined in 
2D square box using numerical simulations based on a 
widely used lattice model of polymers, the BFM. Since 
effects of the competition among configurational entropy, 
bending energy and excluded volume are included, this 
simplified model leads to non-trivial results. By mapping 
the polymer configurations onto the tangent space, we vi- 
sualized changes in conformation of the polymer as the 
bending rigidity was increased from a disordered shape to 
an intermediate elliptical shape and eventually to a spi- 
ral shape (Fig. [3l|4]). We introduced an order parameter 
(Eq. [4|) to characterize the shapes and the distribution of 
the order parameter provided further insight into the na- 
ture of the shape changes. We used a meanfield approach 
to calculate the tangent-tangent correlation function and 
the effective persistence length. 

The combination of numerical and analytical re- 
sults provides a detailed picture of the conformational 
and elastic properties of a polymer under extreme 
confinement. The variation of the persistence length 
with confinement has been studied previously under 
tube confinement where the polymer is unconstrained 
in one dimension. Our results make specific predic- 
tions about the more extreme form of confinement 
where the polymer is constrained in all dimensions. 
The definition of the order parameter, introduced to 
characterize shapes in two dimensions, can be easily 
generalized to higher dimensions and should prove useful 
for characterizing shapes of polymers under confinement. 
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